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I-  ABSTRACT  (Continue  on  roomoo  aide  II  nocooomry  on4  Identity  by  block  mmtbm) 

— ■ — p Investigations  into  the  problem  of  maximizing  array  gain  (DI)  for  situations 
where  ambient  noise  is  the  dominant  interference  have  been  conducted  over  the 
preceding  two  decades  and  possibly  earlier.  These  investigations  have  resulted  in 
mathematical  solutions  having  significant  promise  in  various  practical  applications. 

A closely  related  problem  arises  when  interference  is  dominated  by  reverberation; 
in  this  situation  the  gain  against  reverberation  (the  reverberation  index,  RI)  must  be 
maximized.  The  same  general  mathematical  approach  can  be  applied  to  both  problems-. 


in  that  the  omnidirectional  ambient  noise  limited  condition  can  be  considered  a 
specific  case  of  the  more  general  directional  ambient  noise  limited  condition  and  the 
reverberation  limited  condition. 

Much  of  the  existing  theory  that  has  been  developed  has  been  directed  at  specific 
cases,  thus  leaving  unfinished  a comprehensive  treatment  of  the  general  case.  A 
concomitant  existing  problem  involves  the  determination  of  the  computational  complex- 
ity involved  in  the  mathematical  solution  to  the  described  problems  since  this  complex- 
ity strongly  influences  the  degree  to  which  computational  methods  can  be  applied  to 
practical  problems. 

The  purpose  of  this  paper  is,  therefore,  twofold.  First,  to  solve  the  general 
case,  the  RI  maximization  problem;  second,  to  provide  insofar  as  possible  efficient 
numerical  procedures  for  implementation  of  the  general  solution  in  various  practical 
applications.  These  objectives  are  achieved  by  detailed  mathematical  formulation 
based  on  matrix  theory,  which,  as  will  be  shown,  yields  a complete  procedure  for 
solving  the  general  case.  This  procedure  is  sufficiently  generalized  so  that  the  maxi- 
mization of  DI  is  automatically  included.  Several  illustrative  examples  are  included  to 
demonstrate  the  validity  and  effectiveness  of  this  technique. 


TR  5375 


TABLE  OF  CONTENTS 

Page 


INTRODUCTION  1 

STATEMENT  OF  PROBLEM  2 

MATHEMATICAL  FORMULATION  4 

NUMERICAL  RESULTS  13 

Equally  Spaced  Line  Array  of  Three  Elements 13 

Equally  Spaced  Array  of  Six  Elements  19 

A Cylindrical  Array  of  Fifteen  Elements 36 

FUTURE  RESEARCH  43 

CONCLUSIONS 44 

APPENDIX A-l 

REFERENCES R-l 


i/ii 

REVERSE  BLANK 


TR  5375 


MAXIMIZATION  OF  REVERBERATION  INDEX 


INTRODUCTION 


Within  the  past  two  decades,  problems  with  respect  to  maximizing 
the  directive  gain  of  antenna  arrays  have  been  extensively  studied  and 
some  solutions  have  already  been  produced.  Concerning  the  optimization, 
Uzkov1  presented  an  approach  to  the  problem  of  optimum  directive  antennae 
design  and  obtained  some  important  results;  Tai2  discussed  the  optimum 
directivity  of  uniformly  spaced  broadside  arrays  of  dipoles;  Butler3 
considered  the  optimization  problem  of  nonuniformly  spaced  arrays;  Ma,4 
Cheng,  and  Tseng^’8  presented  an  approach  to  optimize  arbitrary  antenna 
arrays;  Lo,  Lee,  and  Lee7  optimized  directivity  and  signal-to-noise 
ratio  of  an  arbitrary  array;  Gilbert  and  Morgan8  discussed  optimum  de- 
sign of  directive  antenna  array  subject  to  random  variation;  Mermoz9 
showed  that  matched  filtering  techniques  applied  to  particular  condi- 
tions of  crosscorrelation  created  at  the  antenna  by  omnidirectional 
noise  produce  the  optimization  of  the  directivity:  contributions  have 
also  been  made  by  others  such  as  Block,10  Riblet,11  Chu,12  Yaru,13 
Wilmotte14  and  others. 


In  the  sonar  area,  investigations  have  been  conducted  related  to 
maximizing  the  array  gain  of  an  array  dominated  by  ambient  noise. 
Pritchard1 5> 10  discussed  the  maximum  directivity  index  (DI)  of  a linear 
array;  Farran  and  Hill17  studied  wide-band  directivity  of  receiving 
arrays  and  applied  Lagrange  multipliers  to  maximize  directivity  index; 
Edleblute18  introduced  an  optimum-signal-detection  theory  and  used  the 
signal-to-noise  criteria  to  maximize  DI ; Cox19  discussed  the  optimiza- 
tion problems  taking  sensitivity  into  consideration;  optimum  signal 
processing  was  discussed  by  Bryn,20  Nuttall,21  Owsley,22  Lewis,23  and 
others.  These  studies  have  developed  mathematical  solutions  for  the 
optimization  of  DI  with  significant  promise  of  practical  applica- 
tions. A related  important  optimization  problem  arises  when  the  back- 
ground interference  is  nonisotropic,  for  example,  when  it  is  reverbera- 
tion related  or  when  ambient  noise  is  itself  nonisotropic.  Under  such  a 
condition,  the  general  solution,  i.e.,  the  reverberation  index  (RI) 
optimization  is  required.  The  RI  optimization  problem  was  never  con- 
sidered previously.  Adherence  to  the  solution  for  the  general  case 
throughout  all  phases  from  the  initial  problem  formulation  to  final 
modeling  technique  allows  all  aspects  of  an  array's  directivity  to  be 
investigated  using  essentially  the  same  mathematical  procedure.  In 
order  to  encompass  the  already  existing  theory,  some  adaptation  was 
necessary  within  the  work  of  the  general  case  in  the  formulation  of 
mathematical  technique.  As  noted  earlier,  an  important  part  of  the 
development  of  these  techniques  is  related  to  the  investigation  of  the 
degree  of  computational  complexity  necessary  for  the  practical  utiliza- 
tion of  these  techniques.  In  view  of  the  preceding,  the  purpose  of  this 
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paper  is  twofold.  First,  to  solve,  in  general,  the  RI  maximization  pro- 
blem; second,  to  achieve  the  complete  feasible  numerical  treatment  to 
various  practical  applications.  The  entire  problem  is  mathematically 
formulated  and  is  handled  by  means  of  the  treatment  used  for  eigenvalue 
(Gantmacher24)  problems.  This  technique  will  express  the  RI  as  a ratio 
of  two  quadratic  forms  that  determine  an  optimum  set  of  amplitude  and 
phase  excitations  to  yield  the  maximum  RI . This  procedure  is  developed 
with  sufficient  generality  so  that  the  maximization  of  DI  is  included 
automatically.  This  paper  states  the  problem,  presents  a section  of 
theoretical  formulation,  and  includes  a section  of  three  illustrative 
examples.  A later  section  also  discusses  remaining  research  problems. 
Some  conclusions  are  drawn  from  the  theoretical  viewpoint.  The  basic 
computer  programs  that  have  been  developed  to  perform  the  necessary 
computations  on  UNIVAC  1108  computer  are  included  in  the  appendix.  The 
calculation  of  beam  patterns,  RI,  and  DI  is  carried  out  by  the  GBFAM 
(Lee25)  program.  Implementation  of  this  technique  ideally  requires  a 
complete  automation  of  computer  programs.  Such  a package  is  already 
under  our  consideration  and  development,  but  is  not  yet  available  at 
this  reporting  period.  It  is  believed  that  once  this  model  is 
developed,  research  can  proceed  and  meaningful  solutions  can  be 
obtained  to  address  more  complicated,  presently  unanswered  problems. 


STATEMENT  OF  PROBLEM 


Before  we  state  the  problem,  a few  definitions  have  to  be  made. 

We  consider  the  coordinate  system  for  array  structure  according  to 
the  configuration  shown  in  figure  1,  where  6 is  the  azimuthal  deviation 
in  the  horizontal  plane,  4>  is  the  vertical  deviation  from  a horizontal 
reference,  and  Os0£2ir,-Tr/2< <(><*/ 2.  A reference  plane  for  measuring 
relative  phase  at  any  element  of  the  array  is  chosen  to  be  that  plane 
passing  through  the  origin  that  is  parallel  to  all  wavefronts. 

Let  Rj  be  the  amplitude  response  at  the  element,  dj  (0,  $)  be 
the  distance  from  the  jth  element  to  a reference  plane  that  is  taken  to 
be  perpendicular  to  the  direction  (9,$)  of  the  incoming  acoustic  waves, 
and  oj  are  unknown  coefficients  to  be  determined.  Then,  the  beam  output 
function  is  defined  to  be 


n 

V (6 , 4>,  oj)«  l OjRj  exp[i(2ir/X)dj  (0,  ♦,)]. 
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Figure  1.  A 3-Dimensional  Array  Geometry  Structure 


When  aj*l  for  all  j,  the  term  V(6,  $,  1),  conventionally  written  as 
V(6,  $,)  refers  to  an  unshaded  beam  output  function.  V(6,  oj)  is 
used  to  calculate  the  beam  pattern  function,  the  reverberation  index  RI , 
and  the  directivity  index  DI. 

Reverberation  index,  RI,  can  be  expressed  as 


RI*10  logjo 


2ir(sin  b-sin  a)|V(0o,  $0,  aj)|‘ 


17. 


{RV(e,  ♦)>|v(e,  ♦,  aj) | 2cos  * d* 


de 


where  (60,  +Q)  stands  for  steering  direction  angles  and  RV(0,  is  a 
reverberation  function  that  is  used  to  characterize  the  relative  rever 
beration  strength  of  backscattered  energy  from  boundary  or  volume 
scatterers  over  the  angular  intervals  of  interest.  If  we  assign 
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a=-ir/2,  b=ir/2,  and  RV(6,  4>)=1  for  all  0,  <f>,  we  arrive  at  the  expression 
of  directivity  index,  DI,  i.e.. 


4w |V(6  , 0j)|2 

Dl  ■ 10  1 o g j o - 


ri*r*l  2 

JoJ.T'/iln*.  ♦,«.)  |2cos  $ d ♦ de 


Thus  DI  is  mathematically  a special  case  of  what  we  term  RI  and  so  the 
optimization  problem  as  treated  herein  is  equally  related  to  both  RI  and 


DI. 


We  now  pose  the  problem:  Can  aj  be  chosen  such  that  RI  and  DI  are 

at  their  maxima?  This  complete  volume  is  devoted  to  answer  this  ques- 
tion. Details  follow. 


MATHEMATICAL  FORMULATION 


The  beam  output  function,  V(8,  , is  given  by 


n 

V(6,  ♦)-  1 R-iexp  t i(2ir/X)dj  (6,  *)]  , 

j-1  J 


(1) 


where  0 defines  the  azimuthal  deviation  in  the  horizontal  plane,  $ de- 
fines the  vertical  deviation  from  a horizontal  reference,  Rj  is  the 
amplitude  response  of  the  j-th  element,  and  dj  (0,  $,)  is  tne  distance 
from  the  j-th  element  to  the  reference  plane. 

For  a plane  wave  excitation,  the  distance  from  any  element  to  the 
reference  plane  is 

dj(6,  <tO=Xj  cos  a+yj  cos  0+Zj  cos  y,  (2) 

where  (cos  a,  cos  6,  cos  y)  are  direction  cosines  which  define  the 
direction  of  propagation  perpendicular  to  the  incoming  wavefront  and 
(xj,  yj , zj)  stands  for  Cartesian  element  location.  For  a spherical 
wave  excitation, 

♦)-r0-l(Xj-xc)2*(yj-ye)2*(Zj-Z£)2J  l/2.  (3) 
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where  r0  defines  the  distance  of  a point  source  to  the  array  center, 
and 


X£=rQ  cos  <J>  sin  6, 

y^=r0  cos  <(>  cos  6,  (4) 

Z£=-r0  sin  <fi. 

For  details,  the  reader  should  consult  reference  25.  The  directivity 
index,  DI , is  expressed  by  the  formula. 


D I * 1 0 logjgAj'10  log10 


v(e0,  i0) | 


J_f2'r'/2  | v ( e , ♦)  |2cos  <p  d <#>  de 
uwjo  J-ir/2 


(5) 


(Urick26).  The  reverberation  index,  RI,  is  expressed  by  the  formula, 


RI=10  log10A2 


10  log10- 


V(6-,  |2 


(6) 


1 p2*y'b 

2w(sin  b-sin  a)J0  J a{RV(G,  ♦))|V(6,  $)|2cos  ♦ d$  dG 


The  optimization  of  DI  and  RI  is  the  same  as  optimizing  A ^ and  A 2* 

We  shall  attempt  to  optimize  RI  and  DI  in  a systematic  manner.  In 
general  terms,  we  consider  the  reverberation  function,  RV(G,  <t>) , which 
is  used  to  characterize  the  relative  reverberation  strength  of  back- 
scattered  energy  from  boundary  or  volume  scatterers  over  the  angular 
intervals  of  inteiest.  In  the  event,  the  RV(6,  is  a transmitting 
beam  output  function,  the  introduction  of  amplitude  and  phase  shading 
coefficients  to  RV(6,  4>)  should  be  of  interest  for  future  research,  but 
it  is  not  considered  in  our  present  treatment.  Our  approach  is  to 
maximize  A2  and  we  will  show  that  maximization  of  A 2 accommodates  the 
maximization  of  A j. 


Introduce  a set  of  n elements,  {aj},  j = l,  2,  ....  n»  to  V(6,  41). 
Then  we  obtain 


n 

V(e»  ♦.  2 a.R  i exp  |i(2x/ X)  d . ( 6,  ♦))  . 

J J J J 


(7) 
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aj  is  the  j**1  amplitude  excitation  called  shading  coefficients  to  be 
determined.  If  aj=l  for  all  j,  V(6,  $,  1)  = V(9,  $) , which  is  simply 
an  unshaded  beam  output  function.  If  aj  turn  out  to  be  complex,  they 
can  be  put  in  polar  form;  the  magnitude  part  gives  the  amplitude 
excitation  and  the  argument  part  gives  the  phase  excitation.  To  be 
general,  we  consider  the  problem  of  maximizing  A 2-  Once  this  is  done, 
it  is  a simple  matter  to  maximize  A l by  just  setting 


a=-ir/2,  b=ir/2  and  RV(0,  <(>)  = 1, 


which  implies  that  formula  (6)  reduces  to  formula  (5).  Without  loss  of 
generality,  we  drop  the  amplitude  responses  for  V(0,  <|>,  j)  for  formu- 
lation simplicity.  Define  a vector  E t~  be 


r-' 


E={exp[-i(2ir/X)Dj  (0O,  <f>0)]},j=l,  2,  ...  , n,  (8) 


where  Dj(80,  4>0)  are  steering-direction,  (0O,  <ji0) , delays.  Let:  Abe 

a square  matrix  that  is  obtained  by 


A=EE*, 


(9) 


where  * indicates  the  conjugate  transpose,  B be  a square  matrix  whose 
elements  are  to  be  determined,  and  a be  a vector  consisting  of  the 
components  aj . ^2  formula  (6)  can  be  put  in  the  form 


A2 


a *Aa 


a*Ba 


Cio) 


2 

Now  our  aim  is  to  maximize  A2.  Setting  a*Aa  = |V(0O,  <t>0,  aj)  | en- 
ables us  to  determine  the  elements  of  AsCaj^), 

ajk-expli(2ir/X)  (Dj(e0,  ♦O)-Dk(0o,  ♦«,))]  , (11) 
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where  (0O,  $0)  are  steering-direction  angles.  Similarly,  setting 


a* Ba- /*2l,rb{RV(e,*)  } |V(8, ♦,<*.)  |2cos  * d*  de 

2ir(sir  b-sin  ay  0 d a 

gives  the  determination  of  B elements, 

i f'rb 

b j k » — I I exp[i(2w/X)(di-djc)) 

2w(sin  b-sin  a)J o J a (12) 

* { RV  ( 8 , 4>)}cos  $ d$  d6. 

For  a=-ir/2,  b*ir/2,  RV(0,  4>)=1,  (12)  reduces  to 

i r2Y/2 

bik-  — exp  [i (2w/ X)  (d.-d.)]  cos  * d+  de.  (12') 

3 AvJoJ  - */  2 3 k 


We  shall  use  (12')  to  maximize  DI. 

a*Ba  represents  reverberation  power  received  by  the  array;  it  has 
to  be  positive  from  physical  considerations.  This  implies  that  B is 
positive-definite.  Therefore  A 2 can  be  maximized  according  to  the 
following  optimization  theory: 

If  A (a,  a)  and  B(a,a)  are  2 quadratic  forms  in  the  variables  oj 
and  if  B is  positive-definite,  then  the  characteristic  equation  of  the 
form  A(a,  a)  - XB(a,  a)  is 


p(X)=det(A  - XB)=0. 


(13) 


In  view  of  the  matrix  theory  for  eigenvalue  problems  (Gantmacher24) , 
p(X)  has  real  roots  Xj  and  Xj  < Xj+i  for  j=l,  2,  ...  n-1  and  Xj,  Xn 
represent  the  bounds  of  (10),  that  is. 


X1  i 


a*Aa 

a*Ba 


i*n 
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The  right  equality  of  the  above  expression  is  attained  if 


Aa=Xn8a. 


(14) 


Then,  it  is  seen  that 


a *Aa 

X„»A  ■ -E*B_1E>0. 

a*Ba 


(15) 


To  justify  formula  (15)  we  will  prove: 


a*  A a 

(a)  An* for  "optimum"  a, 

a*Ba 

(b)  Xn-E*B_1E>0,  and 

(c)  Xn  is  the  largest  nonzero  eigenvalue  of  (A-XB). 

Proof 


In  view  of  Aa=X  Ba,  we  shall  take  inner  products  of  both  sides 
with  respect  to  a*.  As  a result,  we  obtain  a*Aa  = Xna*Bo,  and  (a) 
results.  Next,  the  conjugate  transposition  of  Aa=XnBa  gives  a* A*= 
Ana*B*.  Since  A,  B are  both  Herroitian,  then  a*Ar=Xna*B.  By  virtue  of 
the  positive-definite  property  of  B,  the  B-l  exists,  then 


a*AB'l-«*Xn. 


Postmultiply  E on  both  sides  of  the  above  equation,  and  in  view  of 
EE*sA,  we  obtain 


(a*E) (E*B“ *E) ■ (o*E) X„. 
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P 


a*E  is  a nonzero  scaler  — ► An  * E*B"l£.  The  positive-definiteness  of 

B— ► 

An-E*B-1E>0. 

(a)  and  (b)  sh  formula  (15)  is  true. 

We  just  provied  that  An  is  a nonzero  eigenvalue  of  (A-AB)a=0. 

Now,  we  want  to  show  that  An  is  the  largest.  Since  E is  an  (n  x 1) 
vector,  EE*=A  must  be  an  (n  x n)  singular  matrix  of  rank  1.  Since  every 
column  of  A is  a constant  times  E (as  in  every  row),  there  is  only  one 
linearly  independent  column  vector  of  A.  Since  B is  positive-definite, 
(A-AB)a=0  can  be  solved  by  B(B“lA-AI)a=0.  A is  of  rank  1,  B-1  is 
positive-definite,  and  (B~lA)  must  have  rank  1.  Therefore,  (A-AB)a=0 
has  only  one  nonzero  eigenvalue  and  the  rest  of  the  eigenvalues  are 
zero  of  multiplicity  (n  - 1). 

Since  An  >0,  An  must  be  the  largest  eigenvalue.  Thus  (c)  is 
f-'  proved  and  justifies  equation  (15).  Making  use  of  (b) , An=E*B_1E,  and 

we  obtain 

AnE  = EE*B*1E-*AnE«AB-1E-*- 
(AB"  1 - AnI)  E»0-*-(A-  A„B)  B“  *E*0  . 

Therefore  B'^E  is  the  eigenvector  associated  with  the  nonzero  An.  In 
view  of  Aa=AnBa,  (A-AnB)a=0,  and  it  follows  that 

oopt"B_1E,  (16) 

which  is  the  required  eigenvector  and  hence  the  optimum  excitation 
vector.  Equations  (15)  and  (16)  constitute  the  complete  solution  to  the 
optimization  problem  and  these  can  be  applied  to  any  sonar  array.  It 
should  be  noted  that  once  B“*  is  found,  it  is  a simple  matter  to 
calculate  RI  and  DI  since  formula  (15)  reduces  to 

a*  Aa 

«E*a. 

a*Ba 

Therefore,  RI  or  DI  can  be  calculated  from 


10  log1()(E*a).  (17) 


T 
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In  order  to  appreciate  the  physical  significance  of  equations  (5)  and  (6), 
consider  the  following  situation:  For  a/a0pt  it  must  be  true  that 

a*  Aa 


a*B  a 

Suppose  we  set  a=E=(exp  l-i(2w/A)Dj ])T,  (j»l,  2,  ...  , n) . This  is  a 
typical  solution  to  the  array  steering  problem  for  the  situation  where 
the  delays  Dj  are  selected  to  provide  an  array  MRA  at  selected  steering 
angles  (0O,  4>0)  . Although  this  solution  does  not  generally  yield  the 
maximum  value  of  DI,  it  possesses  the  more  generally  encountered  charac- 
ter of  array  directivity  computations  and  so  can  be  used  to  relate  the 
term  in  (15)  to  appropriate  physical  quantities.  That  is,  since  A=EE*, 
the  quantity 


a*Aa* (E*E) • (E*E) • (E*E) 2-n2 

is  just  the  square  of  the  number  of  elements  in  the  array.  The  numerator 
of  (15)  is  therefore  the  signal  power  of  the  array  output  when  a per- 
fectly coherent  signal  of  unit  amplitude  is  exciting  the  array.  Also 
for  this  special  type  of  shading  it  is  clear  that  signal  power  is  at 
maximum  (assuming  that  the  aj  are  normalized  such  that  |aj |<_  1 for  all 
j).  Passing  now  to  the  matrix  B we  can  show  that 

sin  | ( 2 w / A ) (dj  -dj,)) 

b j |c- * 

(2«/A)(dj-dk) 

which  is  recognized  as  the  crosscorrelation  between  array  elements  in 
the  presence  of  isotropic  noise.  Thus  B is  the  crosscorrelation  matrix 
of  the  noise.  With  the  choice  of  a previously  made,  a=E— ♦B_1=I  it 
follows  that 


a*Ba-E*E-n . 


Thus  (15)  gives 

2 

DI=10  log10  n'/n=10  log10n, 

and  so  yields  a well  known  result  for  an  array  of  n elements  in  an 
isotropic  noise  field. 

The  central  aim  of  this  investigation  that  exhibits  a departure 
from  the  usual  array  design  procedure,  is  the  maximization  of  the  ratio 
expressed  in  equation  (15)  so  as  to  obtain  best  system  performance. 

Since  a must  differ  from  the  vector  E,  which  we  shall  assign  the  special 
notation  oq  for  reasons  that  will  shortly  become  apparent,  it  can  be 
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seen  that  signal  sensitivity  must  suffer  accordingly  because  the  numer- 
ator of  (15)  will  be  driven  away  from  its  maximum  value.  It  is  a 
desirable  objective,  however,  since  the  result  of  optimization  of  the 
DI  will  usually  justify  the  disadvantage  of  a reduced  array  sensitivity. 
One  must  consider  the  overall  system  performance  to  fully  appreciate 
this  point.  It  may  have  become  apparent  to  the  reader  that  the  term  DI 
as  used  in  this  report  is  more  accurately  termed  array  gain  since  the 
method  when  applied  is  by  no  means  restricted  either  to  an  isotropic 
noise  field  or  to  a perfectly  coherent  signal  excitation.  In  this 
general  usage,  we  shall  define  the  matrix  A as  the  cross-correlation 
matrix  of  the  signal  with  elements  and  the  matrix  B as  the  cross- 
correlation  of  noise  (or  any  interfacing  array  excitation)  with  elements 
bjk.  The  expression  (15)  is  therfore  a maximization  of  array  gain 

a*Aa  Ljrkajak  a^k 
xn« . 

o*Bo  EjEk<ijak  bjk 

Clearly,  if  the  array  is  to  be  optimized  for  an  actual  real-world  situa- 
tion it  is  necessary  to  obtain  sufficiently  accurate  estimates  of  the 
array  crosscorrelations  for  both  signals  and  noise  that  accurately 
depict  the  performance  of  the  system  under  study  and  the  environment 
encountered. 

Further  consideration  of  the  matrix  B yield  some  measure  of  how 
nearly  optimum  the  array  is,  when  in  the  unshaded  state.  In  the  case 
of  DI  (and  also  for  RI)  we  see  that  B can  be  separated  into  the  identity 
matrix  and  a matrix  Bj,  whose  elements  are  (for  all  practical  cases  of 
interest)  less  than  unity  in  absolute  value. 

B«(l*Bj) . 

In  some  practical  cases,  the  spectral  radius  of  Bj  is  less  than  unity. 
This  property  of  Bj  is  sufficient  to  allow  inversion  of  B by  a Neumann 
series  expansion,  that  is. 


B"  '-I-Bj^Bf - 


■ I (-1) 
£ ■ 0 


Thus  using  E*ag, 
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The  series  representation  of  B"*  shows  that  Oopt  can  be  parti- 
tioned into  a convergent  sequence  of  component  terms^  The  first  term 
is  equal  to  oq  and  this  part  of  oopt  as  previously  indicated  is  simply 
the  usual  array-steering  vector.  The  degree  of  optimization  attainable 
beyond  this  base  performance  depends  on  the  nature  of  the  matrix  B^. 

If  Bj=0  no  further  optimization  is  possible  (as  will  be  shown  in  example 
2 in  a later  section)  and  aopt=E,  (the  unshaded  array)  is  already  opti- 
mum. Conversely,  if  Bj  has  elements  sufficiently  removed  from  zero, 
there  is  an  opportunity  for  further  optimization.  This  is  illustrated 
in  the  next  section  where  a superdirective  condition  is  imposed  on  a 
line  array  with  steering  to  endfire  and  element  spacing  <0s)A.  This 
results  in  a large  number  of  off-diagonal  elements  in  the  matrix  B which 
gives  rise  to  meaningful  terms  in  addition  to  the  initial  term.  If  the 
array  element  separation  is  taken  as  an  interger  multiple  of  (^A, 
however,  the  matrix  B.  vanishes  and  the  array  is  optimum  in  the  unshaded 
state.  The  argument  nas  shown  (at  least  in  the  case  of  perfectly 
coherent  signal)  that  it  is  the  nature  of  the  noise  statistics  that 
dictates  the  degree  to  which  array  shading  can  improve  performance. 

Clearly  the  larger  the  crosscorrelation  between  array  elements, 
the  better  the  change  for  improvement.  Viewed  in  another  way,  the 
technique  identifies,  by  purely  mathematical  means,  the  array  shading 
required  to  minimize  the  noise-power  output  of  the  array  while  influenc- 
ing the  signal -power  output  to  some  lesser  degree.  Only  under  some 
very  special  situations,  the  series  expansion  of  B-*  can  be  used 
effectively  for  obtaining  a sufficiently  close  approximation  of  aopt 
without  undertaking  the  horrendous  task  of  inverting  large  matrices  as 
is  the  case  for  large  arrays.  Letting  <Xq=E,  succeeding  estimates  of 
o0pt  can  be  given  by  the  iteration 

av“ao*Blav_i  ''l1  * 

and  the  sequence  can  be  continued  until 


I lav"av-l I l<e»  c arbitrary, 

so  as  to  approach  sufficiently  close  to  the  limit  aopJt.  In  general,  due 
to  the  number  of  arithmetic  operations  involved  and  trie  large  memory 
storage  requirement,  the  series  expansion  is  not  recommended. 
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NUMERICAL  RESULTS 


In  the  sections  below  three  different  arrays  that  illustrate  this 
technique  are  described  and  discussed.  The  calculations  for  the  first 
example  are  performed  by  hand.  Other  computations  are  carried  out  by 
the  UNIVAC  1108  computer  using  GBF.AM  (Lee^)  and  other  programs  that 
have  been  developed  to  implement  this  technique.  In  all  computations, 
single-precision  arithmetic  is  used  and  the  individual  elements' 
response  is  absent  of  consideration. 


EQUALLY  SPACED  LINE  ARRAY  OF  THREE  ELEMENTS 


Array  Description 


On  the  first  example,  we  present  an  endfire  line  array  of  three 
equally  spaced  elements  placed  on  the  z-axis. 


z 


TR  5375 


Input  Considerations 

The  input  parameters  used  are: 
z-spacing=X/4  where  X=wavelength 

Element  locations:  (0,  0,  (-1/2JX),  (0,  0,  (-1/4)X),  (0,  0,  0) 

Frequency  = 4900  Hz 

Plane  wavefront  is  considered 

Steering  direction:  (90°,  90°) 

Shading  Coefficients 

After  normalization  and  separation  of  amplitude  and  phase,  the 
shading  coefficients  of  this  example  are  determined  to  be: 

Amplitude  Coefficients:  (0.66,  1.0,  0.66) 


ui  1 
§ 

i. 


- 
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Phase  Coefficients  (rad):  (2.85,  0,  2.85) 


How  to  obtain  these  coefficients  is  described  next. 


Analytical  and  Numerical  Calculations  of  Maximum  Directivity  Index 


From  our  formulation, 

cos  a = 0, 
cos  6=0, 
cos  y = -zj  sin  4> , 


then  d-;=-Zj  sin  $ for  j =1 , 2,  3.  Using  formula  (12')  to  calculate  bjfc, 
we  find 

1 PVW/2 

b j 1»b22*b  33“—  I I cos  * <!♦  de-1 . 

4wJ  o J.w/  2 

For  jj<k. 


l p»r»/ 2 

b , t ■ — I I exp  |i(2w/X) (d4-d  )l  cos  * d*  de. 
3 4xJo  J -*/  2 3 * 
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Calculations  give 


b l 2*^2  1*1*2  3*b  32*2  / * 


and 


b 


13“b31*°- 


Therefore,  the  B matrix  is  found  to  be 


(1  2 h 

2/ir  1 
0 2/ir 


0 

2/ n 

1 


and 


/ 4 - IT 2 2w 

-4 

1 i 

f 

B-l 

2 TT  -W2 

2 IT 

8 — IT  2 \ 

^ -4  2. 

4 -w 

Since  this  is  an  endfire  array, 

E«  ( exp  I- i ( 2*/  A)  0]  ,exp  I- i ( 2w/ A)  ( - A/4)],  exp  [ - i ( 2w/  A)  ( - A/ 2)  ] ) T 

■(1.  -i,  -1)T- 


And  the  . matrix  is  found  to  be 


* 
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Then, 


opt 


i-  1 


E« 


8-w2 


( 


(8-w2) -2wi 


♦ w2  i 


-(8-w2)  -2wi 


Expressing  aopt  P°lar  form  and  normalizing  it,  we  obtain  ampli- 
tude coefficients  (0.66,  1,  0.66)  and  phase  coefficients  (2.85,  0, 
2.85) . 


/ 16  - 3w  2 \ 2 

lSptAaopt“lV(eo.  ♦o.  «j)  l2*(  8 it2—  ) ; 


aoptBoopt“ 


8 - w ‘ 


Therefore, 


max 


^opt^opt 

aoptBaopt 


7.26 


Max  (DI)-10  log . nD  >8.61  dB 
MO  max 

<n»)un.h.J.d  -4-77  dB 
Improvement  >3.84  dB 


Placing  these  elements  on  different  axes  should  not  alter  the  DI . 
numerical  results  are  shown  by  the  following  table. 


Results 

DI  (dB) 

Unshaded 

4.77 

Analytic  calculation 

8.61 

X-axis 

8.62 

Y-axis 

8.65 

Z-axis 

8.62 

The 
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Beam  Patterns 


We  examined  whether  or  not  changes  occurred  in  the  beam  patterns. 
In  this  example,  both  the  horizontal  and  vertical  beam  patterns  seem  to 
be  improved  somewhat.  The  plots  of  the  endfire  array  on  the  Y-axis  are 
shown  on  page  20. 


EQUALLY  SPACED  ARRAY  OF  SIX  ELEMENTS 


Array  Description 


This 
the  Y-axis 


is  an  endfire  array  of  six  equally  spaced  elements  placed  on 
that  has  the  following  configuration. 


z 


Input  Considerations 


Frequency  = 4900  Hz 

Plane  wavefront  is  considered 

Steering  direction:  (0,  0) 

Since  the  procedure  has  already  been  described  completely  in  the 
previous  example,  the  E vector,  B matrix,  B-l,  and  A vector  calculations 
are  omitted.  We  only  list  the  normalized  shading  coefficients. 
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Using  this  set  of  coefficients  we  obtain  the  following  results 


Max  DI  = 14.54  dB 

(DI)  unshaded  = 7.78  dB 

Improvement  = 6 . 76  dB 

(b)  (1/2) A spacing 

Element  locations  of  this  spacing  are: 

(Xj.  Yy  2 j ) — (0 , tlx,  0),  j = l,  2,  ...  , 6. 

The  determined  coefficients  are: 

Amplitude  coefficients: 

(1.  1.  1,  1.  1,  1) 


Phase  coefficients: 

(u,  0,  w,  0,  n,  0) 


ELEMENT  NUMBER 
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In  this  example,  no  improvement  of  DI  was  found  due  to  element 
spacing. 

(c)  (3/4) X spacing 

Element  locations  of  this  spacing  are: 

(*i.  Yi>  zJ  = (0.  Hllilx,  0),  j-1,  2,  ...  , 6. 

J J J 4 

The  determined  coefficients  are: 

Amplitude  coefficients: 

(0.991,  1,  0.992,  0.992,  1,  0.991) 


a a a 

1 2 3 

4 5 6 

Ul 

O 

3 


ELEMENT  NUMBER 


Phase  coefficients: 

(-.132,  0,  -4.73,  -3.29,  -1.73,  -.412) 


ELEMENT  NUMBER 
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Using  this  set  of  coefficients  leads  to  a small  improvement  of  DI,  0.18 
dB. 


Max  DI  = 7.96  dB 

(DI)  unshaded  = 7.78  dB 


Improvement  = 0.18  dB 


A Comparison  with  Tschebyscheff  Technique:  Using  (1/4) X spacing 
produced  a big  improvement  in  DI  by  this  technique.  This  result  drew 
our  interest  to  an  investigation  if  the  bandwidth  over  which  gain  could 
be  achieved  and  to  a comparison  of  results  against  existing  well-known 
techniques.  The  Tschebyscheff  technique  has  been  selected  for  this 
comparison.  In  comparing  the  DI  improvements  and  beam  patterns,  it  was 
found  that  the  Tschebyscheff  technique  is  inferior  to  this  technique, 
as  shown  in  the  following  table. 


DIRECTIVITY  INDEX: 

(a)  (1/4) X spacing 


Shading 

Function 

Amplitude 

Shading 

Coefficients 

DI 

(dB) 

Improvement 

(dB) 

Gain 

Loss 

Unshaded 

(1,  1,  1,  1,  1,  1) 

7.78 

Tschebyscheff* 

(.3,  .69,  1,  1,  .69,  .3) 

7.04 

0.74 

This  technique 

(.21,  .63,  1,  1,  .63,  .21) 

14.54 

6.76 

*The  Tschebyscheff  coefficients  are  obtained  subject  to  the 
requirement  that  the  side-lobe  level  is  30  dB. 
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(b)  C 1/2) X spacing 


Shading 

Function 

Amplitude 

Shading 

Coefficients 

DI 

(dB) 

Improvement 

(dB) 

Gain 

Loss 

Unshaded 

(1,  1,  1.  1.  1,  1) 

7.78 

Tschebyscheff 

(.3,  .69,  1,  1,  .69,  .3) 

7.04 

0.74 

This  technique 

(1,  1.  1,  1.  1.  1) 

7.78 

0 

(c)  (3/4) X spacing 


Shading 

Function 

Amplitude 

Shading 

Coefficients 

DI 

(dB) 

Improvement 

(dB) 

Unshaded 

(1,  1,  1.  1.  1,  1) 

7.78 

This  technique 

(.991,  1,  .992,  .992,  1, 
.991) 

7.96 

0.18 

Examination  of  Bandwidth 


In  carrying  out  the  complete  computations,  we  used  a fixed  fre- 
quency of  4900  Hz  to  determine  a set  of  shading  coefficients  to  maxi- 
mize DI . The  same  set  of  coefficients  were  also  used  to  calculate  DI 
for  different  frequencies  in  order  to  determine  the  bandwidth.  For  an 
element  spacing  of  (1/4)X,  based  on  a 4900  Hz  center  frequency,  fre- 
quencies ranging  from  3900  to  8000  Hz  were  evaluated  in  100  Hz  incre- 
ments to  examine  the  array  bandwidth.  Associated  horizontal  beam 
patterns  at  f = 4300,  4900,  5800,  7100,  and  8100  Hz  were  selected  for 
presentation.  Fur  a spacing  of  ( 3/4 ) X , based  on  a 4900  Hz  center  fre- 
quency, frequencies  from  4400  to  5400  Hz  were  used;  associated 
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Horizontal  Beam  Patterns 
(a)  X/4  spacing 


« 

S 

? 

a 

2> 

1 

« 


0 0 60  0 12C 

ANGLE  KMgi 

HORIZONTAL  BEAM  PATTERN 
PHI  • 0 0 099 

INDFIRE  LINE  ARRAY  (•)  - X/4  SPACING  - UNSHADED 
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horizontal  beam  patterns  are  presented  at  frequencies  of  4700‘,  4900,  and 
5100  Hz.  A "Frequency  vs  DI"  graph  based  on  the  two  spacings,  (1/4)A 
and  (3/4)A,  considered  above  is  presented.  These  graphs  are  labeled 
figures  2 and  3. 


Discussion 


Let  D(n,  s)  stand  for  the  directivity  of  an  equally- spaced  line 
array  of  n elements  with  element  spacing  s.  Uzkov1  proved  the  following 
important  results  by  means  of  the  orthogonality  properties  of  inner 
products: 


D(n,  s) 


= n when  s = integer  multiple  of  C 1/2) A 
n^  as  s » 0 


(A) 

(B) 


Uzkov's  proof  indicates  that  for  small  element  spacing,  high  directivity 
can  be  expected;  the  technique  just  enunciated  demonstrated  that  such  a 
high  directivity  index  can  be  achieved.  The  results  in  the  previous 
section,  obtained  by  using  a (1/4) A spacing,  clearly  support  the  above 
statements . 

When  a (1/2) A spacing  is  used  in  the  example,  no  DI  gain  is  seen. 
In  such  case,  a DI  gain  should  never  have  been  expected  according  to 
Uzkov's  results.  Mathematically,  under  such  conditions,  our  technique 
determines  matrix  A to  have  unit  elements  and  the  matrix  B to  be  an 
identity  matrix.  Then,  the  largest  eigenvalue  of  B_1  A is  Pn,  which 
equals  n;  this  shows  the  beauty  of  this  technique  and  justifies  Uzkov's 
proof  of  formula  (A) . 

The  maximum  directivity  index  obtained  by  this  technique  by  using 
an  element  spacing  >(1/2)A  appears  to  be  similar  to  those  obtained  by 
Pritchard 15, 16  in  the  process  of  optimizing  the  DI  of  equally-spaced 
line  arrays  of  three,  five,  and  seven  elements.  The  terminology, 
"superdirectivity,"  has  been  used  by  a few  authors.  In  fact,  for  other 
than  half-wavelength  spacing,  this  technique  produces  a narrower  main 
lobe  and  produces  phase  shadings  approximately  the  reverse  of  adjacent 
elements.  These  results  seem  to  match  the  descriptions  given  by 
Urick^6  for  "superdirectivity."  It  is  evident,  in  the  absence  of 
round-off  errors,  our  technique  produces  maximum  DI  (MDI);  this  would 
mean  that  our  MDI  is  also  "super"  in  the  sense  of  superdirectivity. 
Pritchard  concluded  that  for  element  spacings  >(1/2)A  there  is  little 
difference  between  the  maximum  DI  and  the  unshaded  DI . His  conclusion 
confirms  our  findings  for  (3/4) A spacing. 
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Horizontal  Beam  Patterns  (i/4  spacing) 
Tschebycheff 
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Horizontal  Beam  Patterns  (A/4  spacing) 
This  Technique 
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Horizontal  Beam  Patterns  (3X/4  spacing) 


Unshaded 
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horizontal  Beam  Patterns  (3X/4  spacing) 
This  Technique 
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On  the  examination  of  bandwidth,  using  (1/4) A spacing,  we  found  a 
wide  frequency  band  of  4000  Hz;  therefore,  these  results  should  be  of 
interest  to  array  designers.  On  the  other  hand,  with  (3/4) A spacing,  a 
‘.mall  gain  of  DI  is  indicated;  one  should  not  expect  a wide  frequency 
band  as  implied  by  Pritchard. 

To  design  an  array  that  operates  in  the  superdirective  mode,  it  is 
necessary  to  employ  a subtractive  process  on  the  output  of  closely 
spaced  hydrophones.  This  is  often  accomplished  by  reversal  of  the  phase 
of  alternate  hydrophones  in  a line  of  hydrophones  spaced  at  a 1/4  wave- 
length interval;  indeed,  this  technique  does  this.  Unfortunately  the 
subtractive  process,  while  producing  substantial  gains  in  directivity, 
significantly  reduces  the  net  sensitivity  of  the  super-directive  array. 
Moderate  increases  in  directivity  may  be  achievable,  however,  with  a 
moderate  sacrifice  in  sensitivity.  This  was  also  pointed  out  by 
Pritchard. 1 5> 16 

Bandwidth  is  another  limiting  factor  with  superdirective  arrays. 
The  more  superdirectivity  introduced,  the  narrower  the  bandwidth  of  the 
array  with  respect  to  the  bandwidth  that  would  be  available  with  conven- 
tional wider  interhydrophone  spacing. 


A CYLINDRICAL  ARRAY  OF  15  ELEMENTS 


Array  Description 


The  configuration  below  shows  a section  of  a cylinder  that  is 
used  by  this  example.  Array  elements  are  placed  on  three  parallel  row.*, 
each  with  five  elements.  The  separation  angle  is  5°;  the  cylinder 
height  is  1.44A;  the  vertical  spacing  is  0.72A;  the  circumferential 
spacing  is  (n/9)A;  and  the  radius  is  4 ft. 
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Element  locations:  (x^ , y.,  z^)  = (Y  cos  0 j , y sin  0j,  z^) 

0 j“5°  (j-1),  j-1,  2,  ...  , 5 

Frequency:  3500  Hz 

Plane  wavefront  is  considered 

Steering  direction:  (80°,  0°) 

i 1 -1°<*<1° 

RV(0,  ♦)  = 

0 elsewhere 
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Using  this  set  of  coefficients  to  maximize  RI  and  DI  leads  to  the  follow- 
ing improvements. 


Shading 

DI 

RI 

Functions 

(dB) 

(dB) 

Unshaded 

10.29 

5.91 

This  technique 

11.81 

7.19 

Improvement 

1.52 

1.28 

Note  that  any  improvement  is  worthwhile.  In  some  instances  the 
improvement  may  be  greater.  Above  is  only  an  example. 


Discussions 


A noticeable  RI  and  DI  gain  is  produced  as  a result  of  this 
example.  In  carrying  out  the  computations,  we  observed  that  the 
reverberation  function,  RV(e,  <J>) , is  defined  in  a small  sector  that  is 
a part  of  the  main  lobe.  Even  though  the  main  lobe  in  that  small 
sector  is  narrowed  after  the  maximization  is  done,  no  control  can  be 
assured  for  the  beam  patterns  elsewhere  outside  the  sector.  Often, 
when  the  vertical  main  lobe  is  wider  than  the  reverberation  sector, 
undesirable  beam  patterns  can  be  produced.  Under  such  condition,  even 
though  the  RI  is  maximized,  the  DI  may  be  deteriorated.  This  occur- 
rence should  be  considered  unacceptable.  To  remedy  this  situation,  a 
procedure  was  modified  to  improve  the  RI  subject  to  the  requirement 
that  not  only  the  RI  and  DI  can  be  improved  simultaneously  but  also  the 
beam  patterns  can  be  improved.  This  approach  is  tried  to  vary  the 
reverberation  sector  to  cover  the  main  lobe  and  to  improve  RI  and  DI . 
Six  reverberation  sectors  are  tried  for  the  above  example;  the  rever- 
beration function  maintains  the  original  definition.  These  beam 
patterns  are  displayed  below.  Horizontal  beams  all  show  acceptable 
patterns,  but  the  vertical  beam  patterns  deteriorate  with  narrowing 
reverberation  sectors  (indicated  by  [,]). 

The  directivity  index  values  are  also  examined;  they  deteriorate 
with  the  narrowing  reverberation  sector. 
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Beam  Patterns 


Horizontal  Beam  Patterns 


Unshaded 


This  Technique 


ANGLE  (<M«) 

HORIZONTAL  SIAM  PATTERN 


ANGLE  KM9I 

HORIZONTAL  BEAM  pattern 


PHI  - OOcHfl 

CYLINDER  116)  - (BO.  01  - OPT  SHADING 


Vertical  Beam  Patterns 


Unshaded 


This  Technique 


ANGLE  Id*) 

VERTICAL  SEAM  PATTERN 


THETA 


rJ 

a 

f 

\ 

r 

V 

T 

\f 

V 

. 

I 

J 

J 

/ \ 

"7 

/ 

/ 

V 

V 
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'■^^Index 

RI 

DI 

Sector 

(dB) 

(dB) 

Unshaded 

5.91 

10.29 

[-90°,  90° ] 

7.19 

11.81 

[-60°,  60° ] 

7.18 

11.41 

[-45°,  45° ] 

7.19 

9.31 

[-30°,  30°  ] 

7.20 

2.55 

1-15°,  150] 

7.23 

* 

[-1°,  1°  ] 

8.34 

* 

•unreasonable  results 

This  example  represents  a small  array;  the  computation  time  does 
not  impose  a problem.  For  large  arrays,  the  reduction  of  computation 
time  should  be  seriously  taken  into  consideration.  The  saving  can  be 
realized  by  a compromise  on  RI  and  DI  gain.  The  above  results  seem  to 
suggest  that  the  sector  1-60°,  60° ] may  be  a good  choice  for  this 
example.  Now,  the  question  arises:  In  general,  what  is  a reasonable 

sector  to  choose  to  maximize  RI  and  at  the  same  time  to  minimize  the 
computations?  It  seems  that  a good  choice  is  the  smallest  interval 
containing  the  vertical  main  lobe.  In  some  cases,  the  vertical  main 
lobe  is  already  narrow;  then  this  choice  can  save  a great  deal  of 
computation  time. 


FUTURE  RESEARCH 


The  complete  optimization  of  RI  and  DI,  in  fact,  involves  three 
domains  of  analysis:  applied  mathematics,  computer  science,  and 

engineering-physics.  The  solution  of  the  pertinent  problems  in  each 
domain  will  allow  hardware  development  applying  the  subject  technique 
to  be  undertaken.  This  paper  solved  only  the  applied  mathematics 
problem.  Before  the  actual  applications  can  be  carried  out  to  solve 
the  engineering-physics  problems,  it  is  necessary  to  solve  the  problem 
in  computer  science.  It  is  evident  that  the  immediate  need  to  carry 
on  the  research  is  a computer  model  with  desirable  options.  The 
practical  consideration  is  an  automation  of  a set  of  existing  computer 
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programs  plus  options  subject  to  machine  memory  capacity.  Extension  of 
this  model  should  handle  arrays  of  any  size  that  will  be  the  most  gen- 
eral model.  Design  of  such  a model  is  systematic.  For  large  arrays, 
when  the  number  of  array  elements  is  large  enough  to  produce  a large  B 
matrix  that  exceeds  the  allowable  machine  memory  capacity,  the  computa- 
tions are  expected  to  be  extremely  long  and  need  to  be  handled  in  a 
special  manner. 

A complete  examination  of  the  optimum  properties  of  existing 
sonar  arrays  should  be  interesting  and  profitable. 

This  technique  can  be  applied  to  produce  valuable  results  that 
can  lead  to  the  future  design  of  optimal  arrays. 


CONCLUSIONS 

An  effective  technique  based  on  matrix  theory  has  been  developed 
to  apply  an  elegant  theory  to  maximize  RI  and  DI.  Evidently,  the 
numerical  results  demonstrated  its  validity.  Some  engineering  and 
physical  considerations  are  not  examined  in  this  paper.  It  is  believed 
that  this  technique  is  general  enough  to  be  applied  extensively  to  any 
type  of  sonar  array  without  restriction. 

Maximizing  RI  and  DI  is  considered  important.  Also  important  is 
the  minimization  of  computations  for  large  arrays.  Concerning  minimiza- 
tion of  computations,  special  treatment  ought  to  be  considered  for 
special  problems.  For  a large  array,  such  as  equally  spaced  line  arrays, 
the  B matrix  to  be  determined  takes  a special  representation,  a 
Toeplitz  matrix  (Zohar^?) . Existing  special  algorithms  are  recommended 
to  invert  this  type  of  matrix  to  gain  computational  efficiency. 

This  technique  not  only  can  be  considered  useful  for  future  array 
designs,  but  also  can  be  used  to  examine  the  optimum  property  of  pre- 
sent existing  sonar  arrays.  If  an  existing  array  has  already  possessed 
optimum  property,  then  no  improvement  can  be  expected.  This  technique 
will  indicate  so.  In  this  case,  the  matrix  A becomes  a matrix  of  unit 
elements  and  B becomes  an  identity  matrix. 

This  technique  depends  upon  element  spacings  as  well  as  steering 
directions.  Numerical  results  of  equally  spaced  line  arrays  with  small 
spacing  show  superdirectivity.  Cheng,5*®  in  maximizing  directive  gain 
for  circular  and  elliptical  arrays,  concluded  that  for  small  element 
separations  a supergain  situation  results  under  optimum  gain  conditions. 
The  trend  of  numerical  results  seems  to  indicate  that  RI  and  DI  can  be 
maximized  where  the  element  separations  are  very  small. 


‘ V 
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APPENDIX 


A set  of  seven  FORTRAN  V computer  programs  have  been  developed  to 
implement  this  technique.  Their  interdependence  is  described  by  the 
diagram  below.  Calculations  are  performed  on  the  UNIVAC  1108,  using 
single-precision  arithmetic.  Compilation  is  done  by  the  EXEC  8 system. 


1 

1 


MAIN 


XYZ 


EVEC 


•PINT 


FN 


REVB 


1 CGJR 


MAIN  commands  the  control 

XYZ  generates  element  locations 

EVEC  calculates  E vector 

PINT  performs  two-dimensional  numerical  integration 
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FN  computes  the  integrand 

REVB  defines  the  reverberation  function 

CGJR  inverts  Hermitian  matrix  B. 

Outputs  are  a set  of  normalized  amplitude  coefficients  and  asso- 
ciated phase  coefficients. 
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MAIN 


ICO 

1C1 

ICO 


■» 


10 

110 

1C2 

2C1 

3 


PARAMETER  M<=3G 

CO.  .PLla  f)(l,K*Nrx)  *i_'»/(NH)  fiJvd.H)  *V(NX) 

oi  luNSiori  xiiun<)  • yh ( »mk>  *z;i(rim  *X(iir)  .yuir) 

OAT  A fir  ONl  * Zf->’0/0 . 141  j02t»b3b * 1 . 0 * n . 0 / 

OATA  O V/'  a\*  ( 0 . 0 * 0 . 0 ) / 

UAT«  Vr.LC  I Y/49CC  .0/ 

OA,  ...LI  jT  /INFljTb/i»HZ»HVT»LOyH»urMH.LO'.,V,LI.AV,MAWYH»  INL)tXfF»  VELCTY# 
i>  I X YZ * L * ''1 * N * XSPACii  * Y SPACE- * ZSP aC£  * Ti-lEM « E'HIM  , IDA/ 

T ii  J = ..R 

i<E,W<3.  If.i\,TS,LHR  = 1000,ENL>=1000 ) 

■•.L=vKlCTY/A 

CAi_l  XYZ  ( a xYZ  , Xii*  YM* 2H, L *M* *'  * X SPACE * Y SPACE # Z SpACE  * XQ*  YO»  ZO  * l-YM ) 

00  SO  i = l»k_Mi) 

wrtiTE  ( 4 » 1 0 0 1 * Ym(1)  ,ZH(I) 

CO.  i i I.woE 

FORMAT  l 10a,  13,3  (4A,Llb.t»)  ) 

C All  l V...C  t TMlM XH , YH , ZH » LMN *wL*EV) 

Fvj>’.i .a r ( / i ox » vector — */) 

..’KITE  (4*1{J1) 

WHITE  (4»ic3)  (tV  ( J)  * Ozl  »t^N) 
r OHi-iAT  ( 1 C ( / 1 X , 6 (Eib . <3 * La  ) ) ) 

00  10  1 — 1 * l t Ai\| 

<jO  10  0=1*lMN 
i F i * , . o)  GO  To  9 

3(1*0)  =Cf'r*LX  ( Ot'ifc.  * ZERO ) 

v*0  iO  iG 
,'i\  i YH=  i 

vAuL  rSwT  i Hi  1..  $ m V i , LCWf  i * L a MH  * LOW  V * L I * V * 1 . At  I Y 1 1 * Xn  * Y1 1 * LM  * 1 , 0 * 0 * L * f )R  * 
i'O 


wiA  . ft  ii  1 

Call  f X:.T  ( HI <2»HV  T» LOWH , L IMM * LOW  V ♦ L If'V » k’ANYH » Xri * Ym f ZH » I , J , 1 , L#  ID R* 

i’j ) 


01  A*  J)=CFPi_X(A#S) 

CO..T  1 . «UE 
00  110  1 - 1 , LMN 
uO  1 1 0 J=i, LMN 

o(*>»I)=Cf-.PLX(RtAL(lj(I»J)  ) *-AIMAG(u(I*0)  ) ) 
..idlt  (4»lCZ) 

FORMAT  ( /10a  * * 1ST  HOW  UK  13  MATRIX — *//) 
WHITE  i4*iQj)  (b  ( 1 * J)  * 0=1 1 NN’ ) 
rUHf. AT  ( 1H1 ) 

WHITE  (4*2Ql) 

00  0 1=1*  Nt'j 

CUiiTIDOE 

V(i)  = u. 0*0.0) 

CAu_  CGOH  (b*NN*NN» Nil iNn# $1000*  JC»V) 


C-;y  lo'  *©  t • t'"C 
JXnl.at  fa  ly  k-  > . ’ ' 
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DO  o 1 = 1 »Ni\ 

b continue 

DO  *99  1 = 1,  UN 
00  *03  J=i,.\IN 
0V(  I)=;»V(I)+H(1  * J)*EV(J) 
iCo  COmTII.oE 

X ( 1 ) =3 jRT ( KcAl (bV ( 1 ) ) * + a+Alf  Ao (HV ( I ) ) **2) 

Y ( 1 > =A  T AN*  ( A I MAG  ( 0 V ( I ) ) , BuAL  ( HV  ( I ) ) ) 

»«nc  (4^10)  l»OVU)#X  U)fY(i) 

1100  FOKi'iAT  (llA,2(Elb.B»2X)  ) 

PPzX ( I ) *COb ( Y ( 1 ) ) 
oo=a( i )*bii»(Y( n ) 
rfrtliE  (4»H00)  Fr>,00 
*00  CONTINUE 

2i0  FOriMAT ( 2X » 1 3 » 5X » * ( ' ,2(Elb.B,2X) , * ) » # 10X» 'MAGNITUDES* #Elb.8» 1 OX , 
i'A.JoLLs*  »E15.8) 

U=0.0 

uO  11  islf.'xN 

IF ( ( X ( 1 ) ) ,LE  • J)  00  Tu  11 
U=AMb(X(I) ) 
li  co, jrii,oE 

1 a , Oh..,ATl/lbX»A0UH*)//bX#  •bHAOIUG  COO  F 1C  ILNTS  AFTLK  NOrtFAL  I A 1 I 0 
"i»//  10l/iXt3(  Lllb  • U » 4 X ) ) ) 

00  A S i = 1 , i iN 

A(l)=X(I)/vJ 
10  CONTINUE 

..KITE  (4*12)  ( X (o)  » J=1  ,NN) 

30  TO  1 

1000  SToF 
ENO 


Cup/  CTv  . ™L  r-  t 

j*srr. . / 


r c <**•»  ncrl 

. .upiodadfcaiv , 
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XYZ 


bu,r%0071.4fc.  XYZ(lXYd»AH  t YH»ZH , L n'*.*  XSPACE»  YSpACL  #ZSr*ACE»  XO » YO  »Z0 

l 4 i 

mi-dub  lOU  Xil(l)  » YM  ( 1 ) #Zh(l) 
r\.'\U  / 0 , 0 X 7*^3  JiJ  9/ 
ui-iiwL*.'-!*.. 

0 

.itcl.  U)  00  TO  d 
lr (AbPACc  0,0)  oO  TO  3 

l aO=Q 

Auir  uOAT  ( L—l  ) *X bF ACE/ 2.  • 0 
u«J  1 1 = 1 » L 
u>j  i jl  J— 1 r M 
uO  c.  a— 1 1 N 
L X,u- 1 /.  ;+l 

ai  i ( iAO) =r LOAT ( 1-1 ) *XbPAcd-XO 
i 1 CU.-iViwUE 
11  Lui  .T  Xi  iOc. 

X Cui.T  XitOc. 

3 *t- ( 1 j^ACu.  , ud.  0,C)  00  TO  5 

i,  t 0=0 

Yo=G.O 
u*j  A J=  1 » ivi 
<J*J  1H  i = X » L 
i J\J  <i4  A = 1»N 

1 lo=l YG+1 

OUYu)  = (u-l)*Yi,PACE-YO 
.. 4 Cv.'i . i xitud 
.4  wuu  f li  iUii 

4 OOl.)  XNOt 

5 xc  it.br «CL  ,ud.  0,C)  ^0  TO  6 

ido=0 

d^zKLO^r  (,j-i ) *ZbpACE/2,0 
uO  7 A= X t W 

J<j  j.  / 1 = x » u 

Ju  c.  I J=  1 1 M 


dn( idO=Fi_0AT l K— l ) *ZbPACd-ZO 
<.1  OOlJl  iVGd 
17  CGi  jT  l.VJd 
7 wu i jT  I NOc 
o kook;* 
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1 C 
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PINT 
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